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Abstract 

Previous work on soft-gluon resummation for direct photon production is extended to include addi- 
tional subleading logarithmic terms through 0{aag) and some representative comparisons are made 
to experimental results from the E-706 and UA-6 Collaborations. The additional terms are small 
in magnitude, indicating good convergence properties to the level of accuracy calculated. The scale 
dependence remains much smaller than that of the next-to-leading-order calculation. 
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I. INTRODUCTION 



Direct photon production is widely recognized as a process that is potentially important 
in determinations of the gluon distribution function. The next-to-leading-order (NLO) cross 

riQ 

section for direct photon production has been given in Refs. [1|, The role of higher-order 
soft-gluon corrections has also been addressed more recently. Threshold resummation studies 

BttBa-.,....™ 

transverse momentum resummation formalism has been given in Ref. [21]. 

In a previous paper Q we presented analytical and numerical results for the next-to-next- 
to-leading-order (NNLO) next-to-next-to-leading-logarithm (NNLL) soft-gluon corrections for 
direct photon production. Here we increase the accuracy of our previous calculation by including 



additional subleading soft corrections. Our approach follows Re. 
and extends previous work on threshold resummation 

say 
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lich in turn is based on 
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II. NNLO SOFT CORRECTIONS 



At lowest order, the parton-parton scattering subprocesses are q{pa) + g{Pb) liP'y) + i{Pj) 
and q{pa) + <i{ph) lip-y) + diPj)- We define the Mandelstam invariants s = {pa + Pfc)^, 
t = {pa — P-yY, and u = {ph — P'yY, which satisfy S4 = s + t + u = 0at threshold. Note 
that the photon transverse momentum is px = (tu/sY^'^. Here we calculate the cross section 
d^a/d^p^ in single-particle-inclusive kinematics in the MS scheme. The soft corrections to 
the cross section appear in the form of plus distributions 

'ln'(s4/pr)' 
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with / < 2n — 1 at nth order in beyond the leading order, while the virtual corrections 
appear in 6(34) terms. 

We begin with the NLO soft and virtual corrections in the MS scheme. A somewhat different 
notation from that used in Ref. p has been adopted here, so the previously calculated terms 
are repeated here, as well. The corrections to the parton-level cross section, a, can be written 
for either subprocess as 
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where fi^ is the renormahzation scale, and the Born terms are given by 
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where is the charge of a quark of type q, and Cp = {N^ ~ 1)/ i'^Nc) with = 3 the number 
of colors. Also cf = ACp - Ca, cf = Cp + 2Ca, 



of 
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J19 
^2 



[2A) 



where /ij? is the factorization scale, Ca = N^, and (3q = (IIC^ — 2nj)/3, with Uf the number 
of quark flavors. We also define for use below Tl' = — (/3o/4) + 2CfIio.{p'^/ s) and = 
-{3/4:)Cf + {Cf + Ca) HPt/s)- Finally we write c{'^' = + t('^' . For qq -fg we have 
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and Tf = CF[3/2 + ln(p2,/s)] ln(p2./s) - (/3o/4) ln(/iys) +c7 where c'f is defined in Eq. (3.11) 
of Ref . . For qg —>■ 'jq we have 
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Cpln(-^)-CAln 



and Tf^ = [3Cf/4 + Cf ln(-t/s) + Ca ln(-M/s)] \n{pl/s) - (/5o/4) ln(/i|/p2.) + c'f where c'f 
is defined in Eq. (3.8) of Ref. M. 

n 

Note that the NLO q coefficients have also been presented in Ref. [5|]. The notation for C3 
and C2 is the same as before, while the notation for splitting ci into and Ti terms for each 
subprocess is new and useful in presenting the NNLO expressions below. 

Next, we turn to the NNLO soft and virtual corrections in the MS scheme. These corrections 
can be written for either channel as 



For the qq — > ■jg process we have 
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+ [cf cf + (cf )2 - C2 {cff - ^ Tf + ^cf In f ^ ) + 2 Cp 
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Here = (67/18 - tv^/Q) - 5nj/9, G = vrV6, and Cs = 1.2020569 ■ ■ -. The function G 
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denotes a set of two-loop contributions and is given by 
,(2) ^ ^ /7. 22. 299 
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We determine in the virtual corrections Rii^'^s only the terms that involve the renormalization 
and factorization scales, denoted as Rt^'^^^'^s and given explicitly by 
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where (3i = 3ACI/3 - 2nf{CF + 5C^/3) and 



-C''f[-^- + 
For the qg 75 process we have 
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The function G^^^ denotes a set of two-loop contributions P| and is given by 
GfJ = Cl (- A + _ ^Cs) + CfCa (^C3 - ^C2 - ^) 



(2.12) 
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Finally, the terms in _R'?g^7g \\yQ^ involve the renormalization and factorization scales, denoted 
as i?^ 95^79^ aj-g given explicitly by 
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FIG. 1: NLO and NNLO results for direct photon productfon in hadronic collisions compared to pp 
data from the UA-6 Collaboration 

n( at = 24.3 GeV. 

For both processes the coefficients of the 1^3(54), 1^2 (54), and Vii^s^) terms were given pre- 
viously in Ref. 0. The additional subleading Vq{s4) and ^(54) terms presented here are new. 



III. NUMERICAL RESULTS 



In order to show the effect of including the new NNLO terms, the same procedure employed 
in Ref. has been used. First^ a complete NLO calculation of the appropriate cross section is 
performed using a program [14] which employs the phase-space slicing technique described in 



Ref. [15|. The original NLO calculation has been extended to include a complete NLO treat- 



ment of the bremsstrahlung contribution. The Set 2 fragmentation functions of 



161 have been 



used along with the CTEQ6M parton distribution functions [17[. In all cases the factorization 
and renormalization scales have been set equal to a common scale fi which has been chosen 
to be proportional to the photon's transverse momentum. Once the NLO results have been 
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FIG. 2: NLO and NNLO results for direct photon productfon in hadronic collisions compared to pp 
data for the rapidity distribution from the UA-6 Collaboration Ql at = 24.3 GeV. 

obtained, the approximate NNLO contributions can be added to them. Several examples are 
discussed below, comparing the NLO and NLO + approximate NNLO results. 



In Fig. nthe NLO 

are compared to data 



dashed curves) and NLO plus approximate NNLO results (solid curves) 



18^ from the UA-6 Collaboration for proton proton interactions. In each 
case the upper (lower) curve corresponds to the scale choice of /i = (2^^). The pattern 
demonstrated previously in Fig. 6 of Ref. [sf is still found to be true, even after the addition 
of the newly calculated terms presented in this work. The scale dependence of the NLO result 
is greatly decreased by the addition of the NNLO terms. Furthermore, one can see that for the 
scale choice of /i = the NNLO scale dependent terms give a negligible contribution. Even 
with the NNLO contributions, the results lie somewhat below the data at the lower values of 
Pt- 

In Fig. |2lthe rapidity dependence is shown for the UA-6 proton proton data. Again, the scale 
dependent NNLO terms give a negligible contribution for the choice n = pt/2 and the overall 
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FIG. 3: NLO and NNLO results for direct photon productfon in hadronic collisions compared to pp 
data from the UA-6 Collaboration 

n( at = 24.3 GeV. 

scale dependence is greatly reduced when the NNLO terms are added. As noted previously, 
the curves lie below the data over the majority of the rapidity range shown. Of course, this 
distribution is dominated by the contributions from the low end of the px range, so this is no 
surprise, given the results shown in Fig. 

In Fig. El the photon prp distribution is shown for the case of pp interactions and compared to 
data from the UA-6 Collaboration . Whereas the pp reaction is dominated by the qg — > 7g 
subprocess, the pp reaction receives additional significant contributions from the qq —>■ 7(7 
subprocess. Nevertheless, a pattern similar to that in the previous two figures is apparent here 
as well. Note, however, that for the case of fi = pr/2, the NNLO contribution is negative, 
further reducing the scale dependence shown in Fig. El As for the pp case, the band formed by 
the theoretical curves lies somewhat below the data at the lower values of pt- 

Finally, in Fig. (Ha similar comparison is made to data from the E-706 Collaboration [igj ]. 
The same behavior seen in the previous figures is evident here, although the theoretical band 
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FIG. 4: NLO and NNLO results for direct photon production in hadronic collisions compared to data 
from the E-706 Collaboration jl^] at piab = 530 GeV/c. 

now lies significantly below the data at the lower end of the range covered by the data. 



IV. CONCLUSIONS 

We have extended the results of Ref. |^] to include additional NNLO subleading logarithms 
resulting from soft gluon corrections for direct photon production. The additional terms are 
numerically small and the results remain qualitatively the same as in the previous analysis. In 
particular, the reduced scale dependence relative to NLO calculations remains. 
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